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Abstract 

The biocontrol strain Pseudomonas sp. Cab57 was isolated from the rhizosphere of shepherd's purse growing in a field in 
Hokkaido by screening the antibiotic producers. The whole genome sequence of this strain was obtained by paired-end and 
whole-genome shotgun sequencing, and the gaps between the contigs were closed using gap-spanning PGR products. The 
P. sp. Cab57 genome is organized into a single circular chromosome with 6,827,892 bp, 63.3% G-i-C content, and 6,186 
predicted protein-coding sequences. Based on 16S rRNA gene analysis and whole genome analysis, strain Cab57 was 
identified as P. protegens. As reported in P. protegens CHAO and Pf-5, four gene clusters {phi, prn, pit, and hcri) encoding the 
typical antibiotic metabolites and the reported genes associated with Gac/Rsm signal transduction pathway of these strains 
are fully conserved in the Cab57 genome. Actually strain Cab57 exhibited typical Gac/Rsm activities and antibiotic 
production, and these activities were enhanced by knocking out the retS gene (for a sensor kinase acting as an antagonist of 
GacS). Two large segments (79 and 115 kb) lacking in the Cab57 genome, as compared with the Pf-5 genome, accounted 
for the majority of the difference (247 kb) between these genomes. One of these segments was the complete rhizoxin 
analog biosynthesis gene cluster (ca. 79 kb) and another one was the 115-kb mobile genomic island. A whole genome 
comparison of those relative strains revealed that each strain has unique gene clusters involved in metabolism such as 
nitrite/nitrate assimilation, which was identified in the Cab57 genome. These findings suggest that P. protegens is a 
ubiquitous bacterium that controls its biocontrol traits while building up strain-specific genomic repertoires for the 
biosynthesis of secondary metabolites and niche adaptation. 
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Introduction 

Many plant commensal strains have been classified into the 
Pseudomonas jluorescens group, which currently includes more than 
fifty named species [1]. Among the strains in the P. Jluorescens 
group, Pseudomonas protegens strains Pf-5 and CHAO (previously 
called P. Jluorescens Pf-5 and CHAO, respectively) have been used as 
model strains in studies on the biosynthesis of several extracellular 
enzymes such as AprA protease, and secondary metabolites such 
as 2,4-diacetylphloroglucinol (DAPG), pyrrolnitrin (Prn), and 
pyoluteorin (Pit), with antibiotic activity in the rhizosphere [2]. 
These exoproducts contribute to plant protection by these strains 
and other root-colonizing Pseudomonas species with biocontrol 
activity. P. protegens CHAO was isolated from the roots of tobacco in 
Swiss soil which was naturally suppressive to black root rot in 
tobacco caused by Thielaviopsis basicola [3]. The fully sequenced 
strain Pf-5, which is very closely related to strain CHAO, was 
isolated from the rhizosphere soil surrounding cotton seedlings in 
Texas, USA [4,5]. The complete genome data of strain Pf-5 are 
available at the Pseudomonas.com website (http:/ /www. 
pseudomonas.com) and provide detailed annotation data of this 
strain. The genome of strain Pf-5, which is 7.07 Mbp, is the largest 



Pseudomonas genome among Pseudomonas spp. (whose genome sizes 
range from 4.17 Mbp) sequenced to date (http://www. 
pseudomonas.com). The complete genome data of strain CHAO 
have recently been deposited and its size is 6.87 Mbp. Such large 
genomes should confer an environmental fitness advantage on 
these strains. 

The expression of biocontrol factors depends on the Gac/Rsm 
signal transduction pathway ([6] and the references cited therein). 
The Gac/Rsm cascade is initiated by the GacS/GacA two- 
component system. The GacS sensor kinase is known to be 
pliosphorylated and activates the cognate GacA response regulator 
by phosphotransfer during slow growth and in the presence of 
unidentified signal molecules. In P. protegens CHAO, activated 
GacA promotes the transcription of non-coding small RNAs 
(sRNAs) termed RsmX, RsmY, and RsmZ. These sRNAs have 
high affinity for the RNA-binding protein RsniA and also for the 
RsmA paralogue RsmE in strain CHAO [7-10]. The RsmA and 
RsmE proteins repress the translation of genes involved in 
secondary metabolism during trophophase. When RsmX/Y/Z 
sRNAs are induced in idiophase, they relieve the translational 
repression of target genes by sequestering the RsmA and RsmE 
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proteins, thereby allowing the synthesis of secondary metabohtes 
[9,11,12]. Thus, mutants defective in the Gac/Rsm signal 
transduction pathway have a reduced abUity to produce such 
biocontrol factors and to suppress plant diseases [11,13]. Two 
additional sensors termed RetS and LadS were shown to provide 
input into the Gac/Rsm pathway in P. Jluorescens and Pseudomonas 
aeruginosa [14,15]. RetS inhibits and LadS activates the activity of 
the Gac/Rsm pathway and, in this way, both sensors affect the 
expression of target genes and biocontrol factors [15]. 

In aspects of the effective and ecological application of 
biocontrol agents, screening for plant beneficial pseudomonads 
has led to advances in biocontrol research. Moreover, the genomic 
sequence of strain Pf-5 has provided new insights into the 
molecular basis underlying the pseudomonad-mediated suppres- 
sion of plant disease. For example, novel natural products 
including the cyclic lipopeptide orfamide A [16], rhizoxin analogs 
[17] and the insect toxin FitD [18] have been discovered through 
genomics-guided approaches. Comparative genome analysis 
among strains in the P. Jluorescens group has also contributed to 
the understanding of other different strains [19]. P. Jluorescens strain 
A5()6 is now commercially available in the USA and is sold as 
BlightBan A5()6 [20], which is primarily used to suppress bacterial 
disease fire blight on pears and apples. The complete genome data 
of strain A506 revealed distinct features that classified this strain 
into subgroup 3 of the P. Jluorescens group, whereas strain Pf-5 was 
classified into subgroup 1 [19]. 

Although P. protegens has been isolated from North America, 
Europe, and Africa [21,22], litde has been reported on the 
characterization of P. protegens in Asia. In the present study, we 
aimed to screen biocontrol pseudomonad strains relative to strains 
Pf-5 and CHAO from fields in Japan, with expectations that their 
genome data should contribute to the study of core and accessory 
genes of the three strains geographically dispersed. We showed 
that one of the isolates producing DAPG exhibited similar 
biocontrol potential against Pythium damping-oflf and root rot in 
the cucumber to that of strain CHAO. We named this isolate 
Cab57, which was fuUy sequenced and identified as P. protegens. 
Cab57 exhibited typical Gac/Rsm activities such as the expression 
of Rsm sRNA RsmZ and antibiotic production. We also showed 
that the genomic comparison and characterization of each strain 
provides new insights into the identity and diversity of this species. 

Results and Discussion 

Isolation of DAPG-producing Pseudomonads from 
Rhizosphere Soil in Japan 

Approximately 2,800 fluorescent pseudomonads were obtained 
from plant roots. Based on PCR analysis of this collection to screen 
the phirt strain, 48 isolates were selected as candidates of DAPG- 
producing isolates and their production activities were confirmed 
by thin layer chromatography (TLC) (data not shown). Of the 48 
strains, five were PCR-positi\'<' for all of the three other antibiotic 
biosynthetic genes (for pyrrolnitrin, pyoluteorin, and hydrogen 
cyanide) that are typically found in P. protegens. One of the five 
strains named Cab57, which was isolated from the rhizosphere of 
shepherd's purse [Capsella hursa-pastoris (L.) Medik.] growing in a 
field in Hokkaido, was selected for further analysis. 

DAPG production levels by strain Cab57 were determined by 
HPLC. The concentration of DAPG was 248±69 nmol per 
optical density at 600 nm (ODgoo) in the strain grown in nutrient 
broth supplemented with 1 % glucose to ODgoo = 4.0 (one ODgoo 
unit corresponds to approximately 1x10^ cells/mL). The concen- 
tration of DAPG in strain CHAO under the same growth 
conditions was 208 ±62 nmol per ODggo, which suggested that 



strain Cab57 is a prominent producer of DAPG, similar to strain 
CHAO. 

Biocontrol Efficacy of Strain Cab57 

To investigate the biocontrol activity of strain Cab57 in a 
natural habitat, we adopted a cucnmhtr-Fythium ultimum patho- 
system, which enabled us to evaluate plant protection efficacy by 
measuring the root and shoot weights. Strain Cab57 was as 
effective a biocontrol agent as strain CHAO, with 95% confidence 
(Table 1). Furthermore, this strain had a similar root colonization 
capacity to that of CHAO, both in the presence and absence of 
Pythium (Table 1). 

The Genome Structure of Cab57 Showed Properties 
Typical of P. protegens 

We conducted whole-genome sequencing to obtain information 
on strain Cab57. The Cab57 genome is organized into a single 
circular chromosome with 6,827,892 bp (Fig. SI) and 63.3% G-l-C 
content. The genome is predicted to encode 6,186 proteins. The 
general genomic features of strain Cab57 as well as those of strains 
CHAO and Pf-5 are listed in Table 2. 

To identify species, the 16S rRNA gene (position 
125095.126636) of sti-ain Cab57 was applied to a BLAST search 
and was found to be 100% identical to that of P. protegens CHAO 
and Pf-5. Strain Cab57 was also characterized using API 20 NE, 
which revealed the same profile as that of strain CHAO with the 
code 0156557. Strain Cab57 was deposited in the MAFF 
Genebank, National Institute of Agrobiological Sciences as 
MAFF212077. 

Recentiy, the JSpecies program has been commonly used to 

compare two genomes for species definition [23]. The whole- 
genome sequences of P. Jluorescens group strains are publicly 
available. Strains of the P. Jluorescens group used in the genomic 
comparison are as follows; P. protegens Pf-5 [5], P. protegens CHAO 
(http://www.pseudomonas.com), P. chhroraphis subsp. aureofa- 
ciens 30-84 [19], P. Jluorescens PfD-1 [24], P. Jluorescens SBW25 [24], 
and P. Jluorescens A506 [19]; covering all the three sub-clades of P. 
Jluorescens group [19]. The average nucleotide identity calculated 
with BLAST algorithm (ANIb) values, which provides a numerical 
and stable species boundary [23], confirmed the adscription of 
strain Cab57 to P. protegens (Table 3); values higher than 96% were 
found for strains of the same species. Together with ANIb, 
tetranucleotide frequency correlation coefficients (TETRA) rein- 
forced the objective boundary for species circumscription when 
TETRA values >0.99 [23]. As shown in Table 3, both strains of P. 
protegens showed values of >0.999, whereas the other strains 
showed the values of <0.988. A BLASTp search of the three 
strains revealed that strain Pf-5 was most abundant in its unique 
coding serjuences (CDSs) and strain Cab57 shared more CDSs 
with strain CHAO than with strain Pf-5 (Fig. S2). Of the 234 genes 
specific to strain CHAO, 64 genes belong to the regions associated 
with mobile genetic c'l(;ments as mentioned later in Table 4 and 
100 genes encode hypothetical proteins of unknown functions. 

Four gene clusters (hen, pit, prn, and phi) encoding the typical 
antibiotic metabolites in P. protegens were conserved in the genomic 
sequence of strain Cab57 (Gene ID 2622-2624, 2829-2845, 
3743-3746, and 5898-5905, respectively, in Tabic SI), which 
supported the results obtained by screening the five strains that 
were PCR-positive for all four antibiotic biosynthetic genes 
described above. In strain CHAO, the production of the major 
extracellular protease AprA, which functions as a biocontrol 
factor, was also shown to be controlled by the Gac/Rsm signal 
transduction pathway [25,26]. Each of ten genomes of the P. 
Jluorescens group, including P. protegens, has a conserved cluster for 
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Table 1. The suppression of Pythium damping-off and root rot in the cucumber by Pseudomonas protegens CHAO and Cab57. 





Bacterial strain 
added* 


Pythium 
added* 


Surviving plants per 
flask (%)** 


Shoot fresh weight per 
flasl< (g)** 


Root fresh weight per 
flask (g)** 


Colonization by P. protegens (log^o CFU/g 
of root) *** 


None 




100 a 


0.75 a 


0.48 a 


ND 


CHAO 




100 a 


0.73 a 


0.50 a 


7.43 ±0.22 


Cab57 




100 a 


0.73 a 


0.48 a 


7.40±0.18 


None 


+ 


6 c 


0.18 d 


0.10 c 


ND 


CHAO 


+ 


50 b 


0.35 c 


0.24 b 


8.43 ±0.07 


Cab57 


+ 


73 b 


0.43 b 


0.30 b 


8.29±0.10 



*P. protegens strains were added at 10^ CFU per g of soil contained within 100-mL flasl<s (30 g of soil per flasl<), after planting three 92-h-old, sterile-grown cucumber 
seedlings per flask. P. ultimum was added as a millet-seed inoculum at 2.5 g/kg of soil before planting. Plants were harvested after 7 days. 

**Data represent the averages of 10 replicates (flasks containing three cucumber plants) per treatment without P. ultimum and 16 replicates per treatment with P. 

ultimum. Means within the same column followed by different letters (a-d) are significantly different (P<0.05} according to Fisher's protected LSD test. 

***The rhizosphere-stable plasmid pME6031 containing a tetracycline-resistant determinant [49] was introduced as a selective marker into the bacterial strains to 

determine their root colonization capacity in soil. ND, not detected. 

doi:1 0.1 371 /journal.pone.0093683.t001 



this protease [19]. The aprA gene cluster is also conserved in strain 
Cab57 (Gene ID 3232-3236 in Table SI). The reported genes 
associated with the Gac/Rsm signal transduction pathway were 
fully conserved in the Cab57 genome with >99% homology to 
those of Pf-5 (Table S2). 

Other typical gene clusters encoding biocontrol factors were 
also explored (Table S3). Gene clusters for pyoverdine, whose 
product has been commonly reported in strain Pf-5 and CHAO, 
were conserved in the Cab57 genome in four different loci (Gene 
ID 2934-2935, 4179-4198, 4266-4276, and 4286-4288) as 
reported in Pf-5 [27]. Gene clusters for Enantio-pyochelin, which 
has been isolated from strain CHAO and its gene cluster is fuUy 
conserved in Pf-5 [28], and for orfamide A, which has been 
identified in strain Pf-5 from mining oi Pseudomonas genomes [16], 
were also conserved in the Cab57 genome (Gene ID 3642-3651 
and 2162-2164, respectively). Phenazine is an antibiotic produced 
by some Pseudomonas strains (excluding P. protegens Pf-5 and CHAO). 
The homology search of the gene cluster over the entire genome 
suggested that the known pathways for the synthesis of phenazine 
may not be present in the Cab57 strain. 

Rsm sRNA Expression in P. protegens Cab57 

To determine whether P. protegens Cab57 conserved the function 
of the Gac/Rsm cascade, we analyzed the expression of the rsmj(^ 
sRNA gene. As DNA sequences of the upstream region of ww^are 
highly conserved between Cab57 and CHAO strains (four 
nucleotides differ in a 200-nucleotide sequence upstream of +1 
site, with absolute conservation in the upstream activating 
sequence (UAS), and in —35, —10, and +1 sites), the transcrip- 
tional fusion rsm^-lacZ containing the promoter region of CHAO 
was assumed to be apphcable to Cab57 by introducing the 
construct into the strain. As shown in Figure 1, the expression of 



this fusion developed in a growth-dependent manner, which is 
consistent with that previously reported for CHAO [7]. We also 
speculated as to whether knocking out the retS gene (for a sensor 
kinase acting as an antagonist of GacS) in strain Cab57 would 
affect the Gac/Rsm cascade. The retS mutant showed an 
enhanced level of expression (Fig. 1), as reported in CHAO [15], 
confirming the conserved activity of the Gac/Rsm cascade. 

A rets Mutant Exhibited Increased Expression aprA and 
Antibiotic Activity 

The rets mutant had been previously reported to cause 
increased antibiotic production in P. protegens CHAO [15]. In 
CHAO, the Gac/Rsm system has been shown to positively control 
the translational expression of aprA, the gene for the major 
exoprotease [29]. As DNA sequences of the upstream region of 
aprA are known to be highly conserved between Cab57 and CHAO 
strains (no nucleotide differs in a 200-nucleotide sequence 
upstream of the start codon), the translational fusion aprA'-'lac^ 
containing the promoter region of CHAO was assumed to be 
apphcable to Cab57 by introducing the construct into the strain. 
Therefore, we followed the expression of flj6r^-7ac^ fusion in P. 
protegens Cab57 and the rei.? mutant (Fig. 1). The expression of this 
fusion again developed in a growth-dependent manner in Cab57 
and its expression level was enhanced in the retS mutant. 

We also investigated whether RetS would affect the antibiotic 
activity of P. protegens Cab57. In the retS mutant, the production of 
antibiotics with activity against Bacillus subtilis was higher than that 
of the wild type (Fig. 2A), as reported in CHAO. We also tested the 
antibiotic activity of this strain and the retS mutant toward the 
phytopathogenic oomycete Pythium ultimum on a PDA plate. The 
growth of P. ultimum was inhibited by the wild type Cab57 strain, 
and was markedly inhibited by the retf mutant (Fig. 2B). Growth of 



Table 2 



:. General genomic features of P. 













G+C 








P. protegens 


Genome size, bp 


CDS no. 


content, % 


rRNA genes no. 


tRNA genes no. 


Accession no. 


Cab57 


6,827,892 


6,186 


63.3 


16 


68 


AP014522 


CHAO 


6,867,980 


6,115 


63.4 


15 


68 


CP003190 


Pf-5 


7,074,893 


6,108 


63.3 


16 


71 


NC004129 



doi:1 0.1 371 /journal.pone.0093683.t002 
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Table 3. Whole genome comparison of strain Cab57 with those of the Pseudomonas fluorescens group. 



strains within the Pseudomonas fluorescens group 

CHAO Pf-5 30-84 PfO-1 SBW25 ASOS 

Average Nucleotide Identity to Cab57 {%) 98.47 98.22 84.29 81.06 80.55 80.32 

TETRA value to Cab57 0.99992 0.99984 0.98795 0.9256 0.94548 0.94512 



doi:l 0.1 371/journal.pone.0093683.t003 

the phytopathogenic fungus Fusarium oxyspomm was also inhibited 
in a similar manner (Fig. 2C). The pigmentation observed in the 
rets mutant on PDA plates may be attributed to the overproduc- 
tion of secondary metabolites. 

Whole Genome Comparison 

Absence of the rhizoxin analogue biosynthesis gene 
cluster in the Cab57 genome. We performed dotplot analysis 




0 0.5 1 1.5 2 2.5 

Cell population density (OD ) 




Cell population density (OD ) 



Figure 1 . Expression of the rsmZc^M and aprA chad genes in P. 
protegens Cab57. The expression of a rsmZcHAo-'acZ fusion carried by 
plVIE6091 (A) and aprAa^/Ko'-'lacZ fusion carried by pME6060 (B) vj/as 
determined in P. protegens Cab57 wild type (WT, closed circles) and the 
refS mutant (ArefS, open triangles). The expression of these hybrid 
proteins was measured according to LacZ activity. Incubation was 
performed in Erienmeyer flasks as described in the Materials and 
Methods. Measurements were conducted in triplicate. Symbols indicate 
averages and the error bars indicate standard deviations. ODgoo nm- 
optical density at 600 nm. 
doi:1 0.1 371/journal.pone.0093683.g001 



to obtain an overview of genome comparisons among the strains of 
P. protegens (Fig. 3). Compared to the Pf-5 genome, whose size is 
7,074,893 bp, two large segments were absent in the Cab57 
genome (Fig. 3A). These two segments spanning about 195 kb 
accounted for the majority of the difference (247 kb) between the 
Cab57 and Pf-5 genomes. These two segments were also absent 
from the CHAO genome, the size of which was 6,867,980 (Fig. 3C). 

One of these segments was the complete rzx gene cluster (ca. 
79 kb. Fig. S3), which has been reported to encode genes for the 
biosynthesis of analogues of the antimitotic macrolide rhizoxin in 
P. protegens Pf-5 [17]. Rhizoxin was firstly isolated from Rhizopus 
chinensis, a fungus causing a disease in rice seedlings, as its virulence 
factor [30]. Rhizoxin was then found to be produced by 
Burkholderia rhizoxinica, the bacterial endosymbiont of the fungus, 
but not by Rhizopus [31]. The rhizoxin biosynthesis gene cluster (rhi 
gene cluster) was reported in B. rhizoxinica [32], and its fuU genome 
sequencing has revealed that this cluster is encoded on the 
chromosome but not in the megaplasmid carried by this strain 
[33]. However, as the rhi gene cluster is flanked by transposase 
genes, its potential as a mobile region of the genome has been 
suggested. Furthermore, a trace showed this cluster was originally 
on one of the megaplasmids of B. rhizoxinica [33]. The high 
conservation of this cluster between two distant species, P. protegens 
Pf-5 and B. rhizoxinica, suggests the evolutionary traces of 
horizontal gene transfer. In the Cab57 genome, the adjacent 
genes of this cluster (Gene ID 3017 to 3027) were fully conserved 
with >95% amino acid sequence homologies. Although whether 
strain Pf-5 acquired this segment, or strain Cab57 discarded this 
segment after its acquisition remains unclear, such difiFerences in 
antibiotic production potential indicate diversity among strains. 

At the upstream region of this deletion, the fit cluster encoding a 
functional insect toxin reported in P. protegens CHAO [18] was 
conserved between Cab57 and Pf-5 (Fig. S3, Table S3). This 
cluster is also conserved in the distant species, Photorhabdus 
luminescens. However, the transmissibUity of this cluster has not 
yet been demonstrated and its evolutionary origin remains 
vmknown. 

Absence of mobile genomic islands PFGI-1 and -2 of Pf-3 
in the Cab57 genome. A 115-kb mobile genomic island 01 of 
Pf-5 (PFGI-1, [5]) was revealed as another large segment absent 
from the Cab57 genome (Fig. 3 and Fig. S4). This island was 
found to be site-specificaUy integrated into one of the two tRNA 
genes flanking both sides of PFGI- 1 as attL and attR [5] , and its 
corresponding locus in Cab57 represents the tRNA'"^^ gene, which 
is predicted to be a target for the insertion of this island (Fig. S4). 
PFGI- 1 has elements of both temperate phages and conjugative 
plasmids showing similarity to the mobile genomic island 
pKLC102 from P. aeruffnosa C [5,34]. The pil cluster that spans 
over 10 kb is contained in this island and should encode mating 
pHi rather than type IV pHi, and PFGI- 1 also carries putative genes 
involved in conjugal DNA transfer [34] . This pil cluster is shared 
by the pathogenicity islands of P. aeruginosa and genomic islands of 
P. syringae. Other notable genes encoded by PFGI-1 such as a 
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WT 



AretS 



B 



o O 




Figure 2. Effects of the retS mutation on antibiotic activity. 

Antibiotic activities of P. protegens strains toward B. subtilis (A), P. 
ultimum (B) and F. oxysporum (C) were evaluated by the size of the 
growth inhibition zone. Antibiotic activities of P. protegens Cab57 wild 
type (WT) and the retS mutant (AretS) were compared. 
doi:1 0.1 371 /journal.pone.0093683.g002 

cluster of cyo genes, which encode components of the cytochrome o 
ubiquinol oxidase complex, have also been predicted to contribute 



to the survival of P. protegens Pf-5; however, their actual function 
remains unclear [34]. Another genomic island of Pf-5 PFGI-2, 
which spans 16.8 kb, was also absent from the Cab57 genome. 
Both of these islands are also absent from the CHAO genome. 

Phage-related elements and transposon. Six prophage 
regions have been reported as Prophage 01 to 06 in the Pf-5 
genome [34] . We identified four prophage regions as Prophage I 
to IV in the Cab57 genome (Table 4, Fig. 4, Fig. S5, Fig. S6 and 
Fig. S7). Each prophage showed homology to one or two of the 
prophages of Pf-5 and was identified occasionally in different 
integrated sites. For example. Prophage IV, a homologue of 
Prophage 03 of Pf-5, was identified in a different location of the 
Cab57 genome where Prophage 05 was identified in the Pf-5 
genome (Fig. S7). Two large prophages (Prophage III and IV) 
were identified in the genomic region of Gene ID 3451 to 3607 
(Fig. 3A, Fig. S6 and Fig. S7). The homologue of Prophage 06 was 
identified in the CHAO genome in the same integrated sites of the 
Pf-5 genome (Fig. S6 and Fig. S8), reflecting the continuity in the 
syntenic line observed by dot plot analysis in the corresponding 
region (Fig. 3C). The homologue of Prophage 03 was not detected 
in the CHAO genome (Fig. S7). 

The prophage integrated in the mutS/cinA region, which has 
been reported in ten genomes of the P. fluorescens group including 
Pf-5 [19], was also conserved in the Cab57 genome in the same 
position at which Prophage 0 1 was identified in the Pf-5 genome, 
and was subsequently named Prophage I. Prophage I consists not 
only of most of the corresponding prophage of Pf-5 Prophage 01 
(PFL_1209 to PFL_1230), which closely resembles their counter- 
parts from the long contractile tail of the serotype-converting 
bacteriophage SfV oi Shigella jkxneri [34], but also the P2-like tail 
assembly region of Prophage 03 of Pf-5 (PFL_2002 to PFL_2017 
of Prophage 03) (Fig. 4). Such a hybrid structure has also been 
found in a homologous region from P. Jluoresens PfO-l, which 
contains both lambda-like and P2-like tail clusters [34]. The 
genome of CHAO also consists of the homologue of Prophage I as 
a hybrid of Prophage 0 1 and 03 in the same position as that of 
Cab57 (Fig. 4). 

Prophage II shows homology to the 3 ' end of Prophage 04 of Pf- 
5 and its homologue of CHAO, and three ORFs in Prophage II 
showed homology to those of Prophage 06 (Fig S5), suggesting the 
chimeric structure of this prophage. The homologue of Prophage 
02 was absent from the Cab57 and CHAO genomes. The 
homologue of Prophage 05 was absent from the Cab57 genome, 
whereas it was present in the CHAO genome (data not shown). 

Unlike the genome of Pf-5, that of Cab57 possesses two IS 
elements in the genomic regions of Gene ID 1981 to 1998 
(Transposon 1) and Gene ID 3883 to 3893 (Transposon 2) with 
putative transposases (Table S4). Transposon 2 also contains a 
bacteriophage host lysis protein, suggesting it is a mixture of 
transposon and prophage. Transposon 2 was identified in the same 
location at which Prophage 06 of Pf-5 was integrated (Fig. S8). 
Both of these regions were conserved in the CHAO genome, 
whereas several genes including the transposase genes were absent 
from each region (Table S4). 

Overall, the integrated sites of prophages were more restricted 
in the Cab57 genome than in the Pf-5 genome. These strain- 
specific segments of the genome were found in limited chromo- 
somal locations, referred to as regions of genomic plasticity. 

Regions unique to the Cab57 genome. Apart from the 
different large genomic islands and prophages found in each 
genome as described above, we listed regions including at least 
four contiguous ORFs missing in the Pf-5 genome in order to 
identify regions unique to Cab57 (as Cluster SI to S 14 in Table 
S5). One prominent difference in the genome of strain Cab57 
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Figure 3. Syntenic dot plot of P. protegens strains. Dot plot comparison of the Cab57 genome versus the Pf-5 genome (A), the Cab57 genome 
versus the CHAO genome (B) and the CHAO genome versus the Pf-5 genome (C). The blue line represents the regions of similarities between the two 
genomes while discontinuities in this syntenic line represent regions of genomic variations at a given locus between the two strains. Positions of 
phage-related elements and genomic island reported in the Pf-5 genome [34] are marked by arrows (P, Prophage; PFGI, Genomic island). The position 
of the biosynthetic gene cluster for rhizoxin analogs in the Pf-5 genome [17] is also indicated by an arrow (rzx). Regions corresponding to Pf-5 
Prophage 03 (P03) and Prophage 06 (P06) in the Cab57 genome are surrounded by a dashed line in panel A and the expanded view is shown in the 
right. 

doi:1 0.1 371 /journal.pone.0093683.g003 



from that of Pf-5 was the presence of a putative nitrite/ nitrate 
assimilation operon (Cluster S2: Fig. 5 and Table S5). In P. 
aeruginosa PAOl, the corresponding locus PA1778 to PA1786 was 
inferred to encode functions required for a nitrite /nitrate 
assimilation operon, due to homologies with genes reported in 
other bacteria, and its function has already been demonstrated 
[35,36]. For example, the deletion mutant of PA1785, the product 
of which is required for transcription of the nitrate /nitrite 
reductase operon, failed to grow on nitrate as the sole nitrogen 
source [36]. In that study, they found a small RNA termed 
nitrogen assimilation leader A (NalA) in the intergenic region of 
PAOl, which is required for die expression of this operon. The 
NalA homologue was identified in the same locus in the Cab57 
genome. The —24 and —12 boxes of the a'^"^ -dependent 
promoter, both of which are indispensable for promoter activity 
upon nitrogen depletion, were conserved in its promoter sequence. 
Unlike PAOl and Salmonella typhimurium, putative palindromic 
NtrC binding motifs, which are involved in the full expression of 
the promoter activity necessary [36,37], were not found. The nalA 
terminator sequence was conserved as a palindromic sequence 
(Fig. 5). This operon was also present in the CHAO genome. 

P. protegens Pf-5 has been reported to be defective in reducing 
nitrates, which was attributed to the absence of this operon [19]. 
Only three of the ten sequenced strains in the P. Jluorescens group 



(strains 06, Q8rl-96, and Q2-87) possess this operon [19]. The 
presence of this operon in strain Cab57 suggests that this is not a 
core component in the same species; therefore, it may provide 
diversity at the species level. 

Other gene clusters include the following; Cluster S3 encoding 
the penicillin binding protein. Cluster SIO encoding formate 
dehydrogenase. Cluster SI 3 encoding the McrBC restriction 
enzyme, and Cluster S14 involved in the O-antigen (Table S5). 
Cluster Sll and Cluster S12 were identified as secondary 
metabolite clusters by the antiSMASH program. From an 
ecological viewpoint, such gene clusters may also contribute to 
niche adaptation by this strain. Cluster Sll (Gene ID 3624 to 
3635) was located near the pch cluster (Gene ID 3642 to 3651, for 
enantio-pyochelin), and we identified the PFL_3485 homologue 
(Gene ID 3640) which was annotated as the ferripyoverdine 
receptor, between these two clusters (Fig. 6). Cluster Sll possessed 
a specific region (Gene ID 3632 to 3635) that was also absent from 
strain CHAO. One of the genes (Gene ID 3632) was annotated as 
the TonB-dependent siderophore receptor, which suggests differ- 
ent potentials for iron acquisition among the strains. 
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Figure 4. Genetic organization of the region surrounding PropKiage I in P. protegens Cab57 and the corresponding regions in P. 

protegens CHM and Pf-5. The common genes in the three strains are in the same color and conserved genes are connected with gray shading, 
whereas genes that were not common in all strains are colored white. The sizes of genes are not to scale. 
doi:1 0.1 371/journal.pone.0093683.g004 



PLOS ONE I www.plosone.org 



7 



April 2014 I Volume 9 | Issue 4 | e93683 



Genomic Diversity of Pseudomonas protegens 



Materials and Methods 

Ethics Statement 

For the isolation of fluorescent pseudomonads, plant roots were 
sampled on private lands in diverse locations (Hokkaido, Ibaraki, 
Kanagawa, Nagano, Aichi and Nagasaki) in Japan after the land 
owners gave permission to conduct sampling on the sites. None of 
these species is an endangered or protected species. 



a. 

m 
a. 
o 



p E £ 



Isolation of DAPG-producing Pseudomonads from 
Rhizosphere Soil and Selection for Whole Genome 
Analysis 

We isolated approximately 2,800 fluorescent pseudomonads 
based on a method described previously [38]. DAPG-producing 
isolates were pre-screened by PGR using the DAPG biosynthesis 
gene /)A/D-specific primers Phl2a and Phl2b (Table S6). phlD^ 
isolates were incubated on dNBYG medium [39] and DAPG 
production was evaluated by TLC [40]. DAPG production was 
also evaluated in strain Gab57 by HPLG [41]. Other antibiotic 
biosynthetic genes were also explored by PGR analysis using 
specific sets of primers; PRND1/PRND2 for pyrrolnitrin, PltBf/ 
PltBr for pyoluteorin, and PM2/PM7-26R for hydrogen cyanide 
(Table S6). Selection by this method revealed five isolates were 
PGR-positive for all four antibiotic biosynthetic genes and one of 
the five strains named Gab57, which was isolated from the 
rhizosphere of shepherd's purse [Capsella bursa-pastoris (L.) Medik.] 
growing in a field in Hokkaido, was selected for whole genome 
analysis. 
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Bacterial Strains and Growth Conditions 

The bacterial strains and plasmids used are listed in Table 5. E. 
coli and P. protegens strains were routinely grown in NYB (2.5% [w/ 
v] nutrient broth, 0.5% [w/v] yeast extract) and Luria-Bertani 
(LB) medium with shaking, or on nutrient agar plates (4% [w/v] 
blood agar base, 0.5% [w/v] yeast extract) amended with the 
following antibiotics when required; ampicUlin, 1 00 |ig/ mL (only 
for E. coli); kanamycin, 25 |ig/mL; or tetracycline, 25 |ig/mL 
(125 |ig/mL for selection of P. Jluorescens). The inoculation 
temperatures were 28°G for P. protegens and 37°G for E. coli. P. 
protegens was grown at 33°G to improve its capacity to accept 
heterologous DNA (e.g., in electro transformation or in matings 
with E. coll). Bacteria were grown in glycerol-casamino acid 
medium (GGM) for other assays [42]. 
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DNA Manipulation 

Small-scale plasmid extraction was performed with a QIAprep 
spin miniprep kit (QIAGEN); large-scale preparations were 
obtained with a QIAGEN plasmid midi kit (QIAGEN). Ghromo- 
somal DNA from P. protegens was prepared with a QIAGEN 
Genomic DNA buffer set and Genomic-tip (QIAGEN). DNA 
fragments were purified from agarose gels with a QIAquick gel 
extraction kit (QIAGEN). Oligonucleotides used are listed in 
Table S6. 



a 



O) CTl O) CTl 



Whole Genome Sequencing and Assembly 

Whole genome sequencing was performed using Roche 
Genome Sequencer FLX Titanium (for 8 kb long paired-end 
sequencing) and FLX+ (for shotgun sequencing) technology 
provided by Operon Biotechnologies K.K. (Tokyo, Japan). A 
shotgun library was prepared and subsequendy sequenced, 
generating 361,203 reads in 211.0 Mbp of sequencing data. An 
8 kb long paired-end library was also prepared and sequenced, 
generating 394,536 reads in 78.5 Mbp of sequencing data. Co- 
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Figure 5. Genetic organization of the nitrite/nitrate assimilation operons in P. aeruginosa PA01 [36,53] and the corresponding 
regions in P. protegens strains Cab57, CHAO, and Pf-5. (A) ORFs encoding different functional classes of proteins are indicated in each color and 
the ORFs conserved among the strains are indicated in the same color. The percentage of amino acid identity is indicated. Other conserved genes are 
colored gray, and strain-specific genes are colored white. (B) Sequence of the intergenic region between Gene ID 1924 and 1925 in P. protegens 
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Cab57. The putative —24 and — 1 2 boxes of the <j^^ -dependent promoter are shown in bold type and the putative transcriptional start site of nalA is 
underlined. The ct^'' consensus sequence is shown above the sequence. The putative nalA terminator sequence is denoted by horizontal arrows. 
doi:1 0.1 371 /journal.pone.0093683.g005 



assembly of the results of both shotgun and paired-end sequencing 
was performed by newbler 2.6, and were fully assembled into 14 
large contigs defined as > 1 kb. These contigs were then integrated 
only one scaffold, the size of which was 6.8 Mbp with several gaps. 



Gaps between the contigs were filled in by sequencing the PGR 
products using Applied Biosystems 3100x1 DNA Analyzer. High- 
fidelity DNA polymerase KOD Plus™ (Toyobo) was used for 
PGR amplification. The prediction of putative coding sequences 



/ / / / 



^/ ^/ ^/ ^V/ <>^/ ^V/ ^V/ ^ ^ 



Pf-5 



Cab57 



/ / //// 



Clusters 11 



Ton-B dependent 
receptor 



^ ^ »^ ^ n"^ o'^ n'^ 




CHAO 



<Jp^ <^ # # # CT^ 



Ikb 



Figure 6. Genetic organization of Cluster S1 1 in P. protegens Cab57 and tKie corresponding regions in P. protegens strains Pf-5 and 

CHAO. Cluster S1 1 ORFs conserved in Cab57 and CHAO genomes are indicated in blue, whereas those specific to the Cab57 genome are indicated in 
red and yellow. ORFs for the pch cluster are in black. The ORF for the ferripyoverdine receptor (fpvA3) is in green, and that for the TonB siderophore 
receptor is in red. The ORFs conserved among the strains are indicated in the same color. This figure was based on sequence data from the 
Pseudomonas.com website (http://www.pseudomonas.com) and from Youard et al [28]. 
doi:1 0.1 371/journal.pone.0093683.g006 
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Table 5. Bacterial strains and plasmids used in this study. 




strain or plasmid 


Description 


Source or 
reference 


Strains 


Bacillus subtilis 


Ml 68 


Wild type 


C. Keel 


Escherichia coll 


DH5a, HB101 


Laboratory strains 


[50] 


Pseudomonas protegens 


CHAO 


Wild type 


[3] 


Cabs? 


Wild type 


This study 


CabS7retS 


AretS 


This study 


Plasmids 


pCR-Blunt ll-TOPO 


Cloning vector, pUC orl; Km^ 


Invitrogen 


pME497 


Mobilizing plasmid, lncP-1, Tra, RepA(Ts); Ap' 


[51] 


pME3087 


Suicide vector, ColEI replicon. Mob; Tc^ 


[52] 


pME3087retS 


pME3087 containing a BamHI/Hindlll 1.3-kb retS region containing a deletion of the refS gene; Tc^ 


This study 


pME6031 


pACYC177-pVSl shuttle vector; Tc' 


[49] 


pME6060 


Translational aprA'-'lacZ fusion; Tc^ 


[25] 


pME6091 


Transcriptional rsmZ-lacZ fusion; Tc"^ 


[7] 


doi:10.1371/journal.pone.0093683.t00S 



and gene annotation was performed using the Microbial Genome 
Annotation Pipeline (http://migap.lifesciencedb.jp/). The circular 
map that illustrates the general genomic feature as well as its 
comparison with other two strains was plotted by using the 
software BRIG [43]. 

Accession Number 

The complete genome sequence of P. protegens Cab57 has been 
deposited in the DDBJ/EMBL/GenBank database under acces- 
sion no. AP014522. 

Species Identification 

To identify species, the 16S rRNA gene of strain Cab57 was 
applied to a BLAST search. The strain was also characterized 
using API 20 NE strips (BioMerieux Vitek, Inc.) incubated at 28°C 
for 48 hours. Whole genome comparisons of strain Cab57 with 
those of the P. Jluorescens group were performed using the JSpecies 
program [23]. This program is commonly used to compare two 
genomes and uses BLAST software [44]. 

Comparative Genome Analysis 

To compare the whole genome of P. protegens Cab5 7 with the P. 
protegens Pf-5 or CHAO strain, the alignment data between each 
two strains' genome sequences were prepared using the LASTZ 
program (Release 1.02.00), and dot plots were then illustrated 
from the ahgnment data by R (version 3.0.1). Genomic regions 
that were unique to each strain were searched for using the Island 
Viewer program [45] and verified manually. Secondary metab- 
oKte production clusters were examined using the antiSMASH 
program [46]. 

Plant Disease Suppression and Root Colonization Assays 

Flasks containing vermiculite were planted with three cucumber 
seedlings each and treated with Pythium ultimum MAFF425494 
and/ or P. protegens. After seven days of incubation, the biocontrol 



activity and root colonization of each strain were estimated as 
described previously [47]. 

Generation of refS-negative Mutant 

Primers are listed in Supplementary Table S6. An in-frame 
deletion in the chromosomal retS gene of P. protegens Cab57 was 
created as follows. Fragments of 620-bp located upstream and 
640-bp located downstream of the retS gene, respectively, were 
amplified by PGR with the primer pairs RetSUF/RetSUR and 
RetSDF/RetSDR, high-fidelity DNA polymerase KOD Plus™ 
(Toyobo), and genomic DNA of P. protegens Cab57 as a template. 
Each two corresponding fragments were annealed and amplified 
as a 1 .3-kb fragment using the primer pair RetSUF/RetSDR. This 
1.3-kb fragment was cloned into pCR-Blunt II-TOPO™ (In- 
vitrogen). The insert obtained was confirmed by sequencing and 
digested with BamHI and Hindlll. After sequencing, these 
fragments were subcloned into pME3087 cleaved at the BamHI 
and Hindlll sites to give pME3087retS. This plasmid was 
mobihzed from E. coli DH5ot to P. protegens Cab57 by triparental 
mating with E. coli HB101/pME497. Excision of the vector via a 
second crossing-over was obtained after the enrichment of 
tetracycline-sensitive cells, generating the retS mutant. 

P-Galactosidase Assays 

|3-Galactosidase activities were quantified by the Miller method 
[48]. P. protegens strains were grown at 30°C in 50-mL flasks 
containing 20 mL of NYB amended with 0.05% Triton X-100 
with shaking at 180 rpm. Triton X-100 was required to prevent 
cell aggregation. 

Detection of Antibiotic Activity 

The antibiotic activities of P. protegens strains were determined 
with B. subtilis Ml 68, Pythium ultimum MAFF425494, or Fusarium 
oxysporum MAFFl 05034 as the reporter. Cultures of P. protegens 
were adjusted to OD^qo = 1.5 for the assay with B. subtilis, and 
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5-|xL samples were spotted onto the GCM plate. After overnight 
incubation at 30°C, cells were killed by UV irradiation on a 
transiUuminator for 5 min. An overlay of B. subtilis revealed 
antibiotic production by growth inhibition zones. Cultures of P. 
protegens were adjusted to ODroo = 1.5 for the assay with P. 
ultimum or F. oxysporum, 20 \iL samples were streaked around the 
edge of the potato dextrose agar (PDA) plate, and an inoculum of 
P. ultimum or F. oxysporum was transferred to the center of the plate. 
The plate was incubated at room temperature until P. ultimum or F. 
oxysporum reached the edge of the plate. 

Supporting Information 

Figure SI Circular genome map of P. protegens Cab57. 

From the inside out, circle 1 represents the mean centered G+C 
content (bars facing outside-above mean, bars facing inside-below 
mean); circle 2 shows the GC skew (G-C)/(G-I-C); 3 (blue) and 4 
(red) represent comparison of gene content with P. protegens strain 
CHAO and Pf-5, respectively. The outermost ring highlights the 
regions of mobile genetic elements (P, Prophage shown in green; 
T, Transposon shown in purple), and gene clusters which are 
absent from Pf-5 genome (S, shown in black) as described in 
Table 4 and, respectively. The inner black circle shows the scale 
line in Mbps. 
(EPS) 

Figure S2 Venn diagram comparing the coding se- 
quence sets of P. protegens strains Cab57, CHAO, and 
Pf-5. The number of orthologous coding sequences (CDSs) shared 
by ail strains is in the center. The numbers in the non-overlapping 
portions of each circle represent the number of CDSs unique to 
each strain. The total number of CDSs within each genome is 
listed below the strain name. Each set of genes was compared to 
the other two sets using BLASTp, and sequence identity cut-olf 
was set at 60% to identify common genes. Note that due to 
multiple hits, gene numbers did not add up to the strain totals. 
(TIF) 

Figure S3 Genetic organization of Jit cluster (for FitD 
toxin) and rzx cluster (for rhizoxin analogs) in P. 
protegens Pf-5 and the corresponding regions in P. 
protegens Cab57 and CHAO. The conserved genes are colored 
gray, and strain-specific genes are colored white. 
(TIF) 

Figure S4 Genetic organization of the region surround- 
ing the genomic island PFGI-1 in P. protegens Pf-5 and 
the corresponding regions in P. protegens Cab57 and 
CHAO. The conserved genes are colored gray, and strain-specific 
genes are colored white. 
(TIF) 

Figure S5 Genetic organization of the region surround- 
ing Prophage U in P. protegens Cab57 and the corre- 
sponding regions in P. protegens Pf-5 and CHAO. 

Homologous genes are connected with gray shading. The genes 
in prophage are colored white and genes outside prophage are 
colored gray. The sizes of genes are not to scale. 
(TIF) 

Figure S6 Genetic organization of the region surround- 
ing Prophage III in P. protegens Cab57 and the 
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